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k> , 1 Abstract 

C^ I Galaxy hierarchical formation theories, numerical simulations, the discovery of the 

Sagittarius Dwarf Elliptical Galaxy (SagDEG) in 1994 and more recent investigations 
suggest that the dark halo of the Milky Way can have a rich phenomenology containing 
non thermalized substructures. In the present preliminary study, we investigate the 
case of the SagDEG (the best known satellite galaxy in the Milky Way crossing the 
solar neighbourhood) analyzing the consequences of its dark matter stream contribu- 
tion to the galactic halo on the basis of the DAMA/Nal annual modulation data. The 
present analysis is restricted to some WIMP candidates and to some of the astrophys- 
ical, nuclear and particle Physics scenarios. Other candidates such as e.g. the light 
bosonic ones, we discussed elsewhere, and other non thermalized substructures are not 
yet addressed here. 
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2 Introduction 

The DAMA/Nal set-up [1, 2, 3, 4] has exploited the model independent annual modu- 
lation signature over seven annual cycles [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12], achieving 
6.3 a C.L. model independent evidence for the presence of a Dark Matter (DM) parti- 
cle component in the galactic halo. Some of the many possible corollary quests for the 
candidate particle have been carried out both on the WIMP class of DM candidate 
particles with various features and increasing exposures and on keV-range pseudoscalar 
and scalar DM candidate particles (to which experimental activities, applying what- 
ever rejection technique of the electromagnetic component of the counting rate, are 
blind). Various possibilities for the candidate and the interactions have also been 
discussed in literature by e.g. [13, 14, 15, 16, 17, 18, 19]. 

Many of the uncertainties and assumptions affecting whatever kind of model de- 
pendent result in the field (such as e.g. corollary quests for candidate, exclusion plots, 
and -in the case of indirect investigation experiments - the determination of evidence 
itself, of the parameters or of the limits, etc.) have been discussed in some details e.g. 
in ref. [3]. 

Here we will make a preliminary investigation on the effect of DM stream contribu- 
tions in the galactic halo restricting the presentation to some solutions for the WIMP 
class already discussed e.g. in refs. [3, 4] and - as regards the stream contribution - 
to the SagDEG case, which has been already addressed in a different way in literature 
e.g. in ref. [20, 21]. Other interesting cases and candidates will be further addressed 
in this light in future works. 

As known, DAMA/Nal exploited the effect of the Earth revolution around the 
Sun on the DM particles' interactions in the target-material of suitable underground 
detectors. As a consequence of its annual revolution, the Earth should be crossed by 
a variable flux of DM particles along the year. In particular, the expected differential 
rate as a function of the recoil energy, dR/dEn (see ref. [3] for detailed discussion), 
depends on time owing to the DM particle velocity distribution in the laboratory 
frame, f{v\vQ{t))] here Vi^it) is the Earth's velocity in the galactic frame as function 
of time. 

This method offers an efficient model independent signature, able to test a large 
interval of cross sections and of halo densities; it is named annual modulation signature 
and was originally suggested in the middle of '80 by [22] . 

In particular the expected counting rate averaged in a given energy interval can be 
expressed by the first order Taylor approximation: S{t) ~ 5o -f SmCosuj{t — t^) with 
the contribution from the highest order terms less than 0.1%. Sq is the unmodulated 
term, S^ is the modulation amplitude, lu— 27r/T with T=l year and ioj time when 
the expected counting rate is maximum, depends on the adopted halo model and on 
possible non thermalized contributions. In particular, for halo models with velocity 
distribution isotropic in the galactic frame io is roughly June 2'"^ when the Earth 
velocity in the galactic frame is at maximum. In the present paper we use the following 
prescription for the velocity distribution in the laboratory frame once fixed the halo 

model: ptot X /(w|we(i)) = Phalo X fhalo{v\V(Bit)) + Pstream X f stream iv\vQ{t)) , whcrC 

the two contributions of the DM particles in the dark halo and the DM particles in the 
stream have been pointed out. Here ptot, Phaio, Pstream are the DM particle densities 



and fhaio and f stream are the velocity distributions of the two components normalized 
to one. 



3 SagDEG phenomenology 

Since the discovery of the SagDEG in 1994 [23], it has been argued that the dark 
halo of the Milky Way can have a rich phenomenology containing non thermalized 
substructures. This hypothesis is also supported by galaxy hierarchical formation 
theories [24] and by some numerical simulations [25]. Additional interest is offered by 
the observation of other satellites of the Milky Way, such us the Canis Major in 2003 
[26], and satellites of other near galaxies, like the stream {v ~ 300 km/s ) discovered 
in our nearest neighbouring "twin" galaxy M31 [27]. In 1998 it was found that the 
SagDEG orbits the Milky Way Galaxy in about 1 Gy, having passed through dense 
central region of our Galaxy at least about 10 times during its life. This has been 
interpreted as an indication of presence of DM that with its gravity has prevented the 
disruption of the SagDEG [28]. 

Suitable DM direct detection experiments can provide interesting information 
about the local halo structure, investigating the presence of non thermalized dark mat- 
ter fluxes, as in case of the tidal stream of a dwarf satellite galaxy passing through the 
solar neighbourhood. SagDEG is the best known satellite galaxy of the Milky Way 
crossing the solar neighbourhood; here the consequences of its dark matter stream 
contribution to the galactic halo are analyzed on the basis of the DAMA/Nal annual 
modulation data [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. In fact, e.g. - as reported in ref. 
[4] - the presence of DM streams in the Galaxy would induce a (slightly) variation 
of the phase value of the modulated component of the signal and its variation with 
energy. Consequently, also the ^ ratio undergoes a change depending on the energy 
window and on the stream properties ^. 

The SDSS and the 2MASS surveys [29, 30] have traced the tidal stream of the 
SagDEG; two streams of stars are being tidally pulled away from its main body and ex- 
tend outward from it. The leading tail can shower matter down through the solar neigh- 
bourhood and considerations, based on the (very uncertain) M/L ratio, suggest the 
allowed density in the SagDEG tail, psgr, to be of the order of (0.001 — 0.07)GeV cm^^ 
corresponding to about (0.3 — 23)% of the halo local density [20]. 

Fundamental information, in order to investigate effects correlated with the pres- 
ence of such tidal streams, is the value of the mean velocity of the stream, its direction 
and its velocity dispersion. Despite the fact that SagDEG is the best known satellite 
galaxy of the Milky Way crossing the solar neighbourhood, these quantities are not yet 
well defined and a large number of related investigations can be found in literature. 
In particular, in this paper we have taken into account both the values of ref. [20] 
(derived from the analysis of eight clump stars - from Chiba and Yoshii catalogue [31] 
- attributed to the SagDEG tail) and the values of ref. [32] (based on a SagDEG 
simulation model). 



^In particular, the -^^ ratio would increase or decrease with the respect to absence of stream in 
the galactic halo depending on the local direction of the stream in the halo. 



In the following, we will use a right-handed reference frame with the x axis towards 
the Galactic Center, the y axis towards the direction of Galaxy rotation and the z axis 
towards the Galactic North Pole. 

To account for the determination of ref. [20] , the SagDEG stream has been modeled 
as a DM flux with mean velocity in galactic coordinates given by: 



Vs* = iV^,Vy,V^)^ (-65 ± 22, 135 ± 12, -249 ± 6) km/s. 



(1) 



Here 1 a error has been reported for each velocity component. In addition, to account 
for the determination of ref. [32] we have also considered the following cases : 



and 



Vsph = {V^,Vy,V,) = (-86 ±14, 69 ±3, -384 ± 1) km/s, 



Vobi = iV,,Vy, V,) = (-57 ± 8, 79 ± 3, -395 ± 1) km/s, 



(2) 



(3) 



for the spherical and the oblate halo models of ref. [32] , respectively. These two last 
stream mean velocities have been derived by considering for each halo model the ~ 
100 configurations nearest to the Sun within a distance < 2.5 kpc (see Fig. 1). It 
is worth to note that the prolate model of ref. [32] has not been considered in the 
following since no configuration is present in the solar neighbourhood. A graphical 
representation of the three stream mean velocity sets, considered in the present paper 
for the SagDEG tidal stream, is shown in Fig. 2. 






Figure 1: SagDEG simulation models for spherical (left), oblate (center) and prolate 
(right) halo potentials; data taken from ref. [32]. In each panel the circle pointed by 
the arrow selects the Earth position and the configurations considered in this paper 
for the evaluation of the used mean velocity values: Vsph and Voti- We note that no 
configuration is present in the solar neighbourhood for the prolate model. 



As regards the velocity dispersion to be associated with each one of the three 
considered stream parametrizations given above, we derive for our reference frame, 
respectively: 



{(7x,<Jy,<Jz)s* = (62,33,17) km/s, 
i(Tx,a-y,<Jz)sph = (60,19,8) km/s 



(4) 
(5) 




Figure 2: Graphical representation of the three stream mean velocity sets studied in 
the present paper for the SagDEG tidal stream. 



and 



{(Tx,cry,crz)obi = (59,23,9) km/s. 



(6) 



The {crx,cFy,cFz)8* is taken from ref. [20], while {ax,cry,az)sph and {ax,cry,az)obi have 
been calculated for each model as r.m.s. values of the about 100 configurations in 
the solar neighbourhood. As it can be observed, for the three considered SagDEG 
stream models, the velocity dispersions are quite different and significantly non- 
isotropic. Notwithstanding, in the following for simplicity the velocity distribution of 
the SagDEG stream in the solar neighbourhood has been approximated by an isotropic 
Maxwellian distribution in the locally comoving frame of SagDEG. The investigation 
of the effect of non-isotropic distributions is not considered in the present paper and 
can be addressed in future. 



3.1 A full example of the SagDEG effect on the annual modu- 
lation signature for a given scenario 

In this subsection we show for template purpose a complete example of the effect 
induced by the presence of a DM stream in the solar neighbourhood considering a 
particular scenario. 

As known, the Earth velocity in galactic coordinates can be expressed as following: 



w©(i) = VLSR + vq + VEarth{ei sinA(i) - 62 cos A(t)) 



(7) 



where vlsr — (0,220 ± 50,0) km/s (the quoted uncertainty is at 90%C.L.) is the 
velocity of the Local Standard of Rest; vq = (10.00, 5.25, 7.17) km/s is the Sun pecuhar 
velocity here taken from ref. [33] and Vsarth is the mean orbital velocity of the Earth 
(~ 29.8 km/s). The ei, 62 versors and the A(i) function are [21]: 

ei = (-0.0670,0.4927,-0.8676), 

62 = (-0.9931,-0.1170,0.01032), 

\{t) ^uj{t- 0.218). 

Here uj = 27r/T with T = 1 y, i is the time in years starting from January 1st and 
0.218 y is the spring equinox (March 21). 

The velocity distribution of the SagDEG DM particles in the laboratory frame can 
be written as ■^: 



fsariv) = ^^ e <^^^ , (8) 

where the mean velocity of the SagDEG DM particles in the laboratory frame is: 
Vsgr,e{t) = Vsgr ~ «©(*); Vggr wiU bc in turn either Vs* or Vsph or Voti- Finally, for 

each Vsgr, the vo,sgr — \ h^^sgr parameter is assumed in the following to be cither 20 
or 40 or 60 km/s. 

The |wsgr,©(i)| reaches its maximum value at time to.sgr defined by: 

CO.9 [A(to,s3r)] = — F=^==, sin [A(io.sgr)] = , ^^ „ , (9) 

with tti — Ci ■ {vsgr — vj^sR ^ Vq). Therefore, the mean velocity of the SagDEG DM 
particles in the laboratory frame would be maximum around January 10*'* — 14*'', 
depending on the considered SagDEG velocity set. We remind that in absence of 
SagDEG contribution (that is, fhaio only contributes to the total DM particles velocity 
distribution), tg is expected to be roughly at the 152.5 day of the year (~ June 2"''). 
Hence, the net effect of a SagDEG tail contribution to the local halo density is a shift 
of few days (towards January) in the expected phase of the signal. For the sake of 
completeness, the ^ ratio is not expected to be enhanched when a SagDEG stream 
is included, due to the nearly opposite phases between the two extreme cases. 

As examples, in Fig. 3 the phase, io, of the modulated component of the signal 
is plotted - for some models with the inclusion of the SagDEG stream and for some 
reference WIMP masses - as a function of the detected energy in Nal(Tl) detectors. 
In particular, in the given examples, simple assumptions have been adopted: 

i) Vsgr = (^65, 135, —249) km/s, that is T^* at its central value; 

ii) wo,sgr = 40 km/s; 

iii) vj^sR — (0, 220, 0) km/s , that is at its central value; 

iv) SagDEG tail DM density psgr — 0.04 x phaio', 

v) Galactic halo model: NEW (a = 1, /3 == 3, 7 = 1, a = 20 kpc) (A5 of [3]); 



^We note that in the following, the quantities related to SagDEG are marked as sgr. 
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Figure 3: Examples of the effect of the SagDEG tail, modeled as given in the text, on 
the expected annual modulation signature in Nal(Tl) detectors. In each panel a plot 
of the phase (to) vs the detected energy (E) is shown for a given WIMP mass and 
for the given assumptions (see text). Dotted line: phaio = 0.74 GeV cm~^ (maximum 
value allowed for the adopted halo model) [3]; dashed line: phaio — 0.33 GeV cm^'' 
(minimum value allowed for the adopted halo model) [3]; solid line: absence of SagDEG 
contribution, that is Iq ^ June 2"''. The effect of a possible SagDEG contribution is 
to slightly shift the phase tg towards lower values at low recoil energy. 



vi) WIMP DM candidate with dominant Spin Independent coupling (cr ex A^); 

vii) Form factors and quenching factors of ^'^Na and ^^^I as in case A of ref. [3]; that 
is, the most cautious Helm form factor, the mean nominal values for the parameters 
of the nuclear form factors and for the measured ^^Na and ^^^I quenching factors are 
assumed. 

For the sake of completeness, we remind that the DAMA/Nal results (107731 
kg-day exposure) provide to — (140±22) day averaged in the (2-6) keV energy window; 
at present level of sensitivity - as it can be seen in Fig. 3 and extensively in the following 
sections - it is consistent both with presence and with absence of SagDEG contribution. 
As discussed in [4] larger exposures, which will be available in near future thanks to 
the presently running DAMA/LIBRA set-up [34], will offer the possibility of more 
stringent constraints. 



3.2 Investigating the effect of a SagDEG contribution for 
WIMP cases 

In order to further investigate the effect of the presence of a SagDEG stream, we wiU 
follow in this section the same approach already exploited in ref. [3, 4], where the 
simplier case without SagDEG contribution was considered. In particular we have 
considered here the WIMP class of candidate particles in the general case of mixed 
SI(Spin Independent)&SD(Spin Dependent) coupling and the two subcases of pure SI 
and pure SD couplings. 

In fact, since the ^"^Na and ^^^I are fully sensitive to both SI and SD interactions 
the most general case is given by a four-dimensional volume {m\Y, C'^'s/? ^'^SDi ^)i 
where miy is the DM particle mass, £, is the ratio between the local density for the 
considered candidate and the local DM density ptot, crsi is the SI WIMP-nucleon cross 
section and aso is the SD WIMP-nucleon cross section according to the definitions 
and scaling laws considered in ref. [3]; tg9 is the ratio between the effective coupling 
strengths to neutron and proton for the SD couplings (6 can vary between and 
tt). In the calculation the same galactic halo models and associated parameters as 
in ref. [3] have been considered as well as the uncertainty on the value of the local 
velocity vo = (220 ± 50) km/s (90% C.L.). For the case of the SagDEG stream 
description we have considered the three possibilities for the velocity (T^*, Vgph and 
Vobi at their central values) and three possible velocity dispersions {vo,sgr — 20, 40, 60 
km/s), obtaining nine different cases for each fixed halo model. Moreover, in this 
section we consider that the SagDEG contribution cannot exceed ~0.1 GeV cm^'^, as 
suggested by M/L ratio considerations of ref. [20] . 

The results presented by DAMA/Nal on the corollary quests for the WIMP 
candidate particles over the seven annual cycles are calculated here and elsewhere 
(e.g. [3, 4] and references therein) taking into account the time and energy be- 
haviours of the single-hit experimental data. For this purpose, the likelihood func- 
tion LI"' {mw, S.CTSI, S.o'sD, 0) is constructed for any fixed SagDEG velocity set and 
velocity dispersion (cumulatively labeled here as is) and for all the considered model 
frameworks (cumulatively labeled here as im , running on the galactic halo models and 
on all the other parameters involved in the calculation). In particular, the likelihood 
function requires the agreement: i) of the expectations for the modulated part of the 
signal with the measured modulated behaviour for each detector and for each energy 
bin; ii) of the expectations for the unmodulated component of the signal with the re- 
spect to the measured differential energy distribution; iii) for WIMP candidate (since 
ref. [8]) also with the bound on recoils obtained by pulse shape discrimination from 
the devoted DAMA/NaI-0 data [35] . The latter one (used when WIMP candidates are 
considered) acts in the likelihood procedure as an experimental upper bound on the 
unmodulated component of the signal and - as a matter of fact - as an experimen- 
tal lower bound on the estimate of the background levels by the maximum likelihood 
procedure. Thus, the C.L.'s, we quote for allowed regions, already account for com- 
patibility with the measured differential energy spectrum and, - for WIMP candidates 
- with the measured upper bound on recoils. In particular, in the following for sim- 
plicity, the results of these corollary quests for the candidate particle are presented in 
terms of allowed regions obtained as superposition of the configurations corresponding 
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Figure 4: Examples of slices of the four-dimensional allowed volume in the {£,(tsi , S.o'sd ) 
plane for some mw and 9 values in the considered scenarios. The shaded regions have 
been determined for no SagDEG contribution [3] , while the areas enclosed by the lines 
are obtained by introducing in the analysis the SagDEG stream with DM density not 
larger than 0.1 GeV cm""^. The nine considered possibilities for the SagDEG stream 
veloctity (V^* (blue), Vsph (black), Vou (red)) and wo,sgr dispersion (20 km/s (dashed), 
40 km/s (solid) and 60 km/s (dotted)) have been reported. 



to likelihood function values distant more than 4ct from the null hypothesis (absence 
of modulation) in each of the several (but still a limited number) of the possible model 
frameworks considered here. Obviously, these results are not exhaustive of the many 
scenarios possible at present level of knowledge (e.g. for some other recent ideas see 
[36, 37]) and larger sensitivities than those reported in the following would be reached 
when including the effect of other existing uncertainties on assumptions and related 
parameters [3, 4]. 

For the general case of a WIMP with mixed SI&SD coupling, one obtains a four- 
dimensional allowed volume^. Since a full picture of this result is not possible in 
the practice. Fig. 4 shows some slices of the four-dimensional allowed volume in the 
plane ^asi vs ^aso for some of the possible mw and 6 values. The filled areas show 
the case without SagDEG contribution and, therefore, have already been reported 
and discussed in ref. [3, 4] (some different slices are also shown here), while the 
areas enclosed by lines show the cumulative effect of the possible SagDEG stream 
contribution in various cases (see figure caption). 
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Figure 5: Region allowed in the {S.(7 s i , fnw) plane in the considered scenarios for 
pure SI coupling. The filled region has been determined for no SagDEG contribution 
[3, 4], while the areas enclosed by lines are obtained by introducing in the analysis the 
SagDEG stream with DM density not larger than 0.1 GeV cm~^. The nine considered 
possibilities for the SagDEG stream veloctity (Vg* (blue), Vsph (black), Vom (red)) and 
Vo^sgr dispersion (20 km/s (dashed), 40 km/s (solid) and 60 km/s (dotted)) have been 
reported. 



The purely SI subcase^ is shown in Fig. 5, while in Fig. 6 some slices of the 
three-dimensional allowed volume {mw, ^cfsd, 6) for the purely SD case are given. 

^Is worth to note that - for example - experiments using either nuclei largely insensitive to SD 
coupling (as e.g. ""'Ge, ""^'Si, "''*Ar, "''*Ca, ""'W, "°*0) or nuclei in principle all sensitive to such 
a coupling but having different unpaired nucleon (neutron in odd spin nuclei, such as ^^^Xe, ^^'^Xe, 
^■^^Te, ^'^Ge, ^^Si, ^^'^W) with the respect to the proton in ^^Na and ^^''l cannot explore most of the 
four-dimensional allowed volume. 

®We remind that no direct comparison is possible also among results on purely SI coupled WIMPs 
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Figure 6: Examples of slices of the three-dimensional allowed volume in the 
{£,o'SD,'mw) plane for some 9 values in the considered scenarios and for pure SD 
coupling. See Fig. 5 for the meaning of the regions. 



The filled areas and the areas enclosed by the lines have the same meaning as before. 

As it can be observed, the inclusion of the SagDEG stream and of the related 
uncertainties significantly modifies the allowed volumes/regions; the role appears larger 
mainly for larger WIMP masses. 

It is worth to note that other streams can potentially play more intriguing roles 
and will be investigated in near future, such as the Canis Major [26]. Moreover, other 
kinds of streams as those e.g. arising from caustic halo models [37] can also play a 
significant role in the corollary investigations for the candidate particle with whatever 
approach and for comparisons. 

This approach to the problem allows to make some cautious constraints on the 
SagDEG stream in the galactic halo on the basis of the measured DAMA/Nal annual 
modulation data. We will discuss some of the implications of the presented results in 
the next section. 



achieved by using diflFerent nuclei, although apparently all the presentations generally refer to cross 
section on nucleon. For some discussions on generalities and comparisons see e.g. [3, 4, 38] 
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4 Constraining the SagDEG stream by DAMA/Nal 

As mentioned, the high exposure of DAMA/Nal can ahow to obtain some prehminary 
information about the presence of substructure in the halo as the SagDEG stream. 
For this purpose, the hkehhood ratio function has been used as statistical analysis 
approach in the investigation on the SagDEG parameters with the respect to all the 
others involved in the calculation. 

In particular, fixing: i) a SagDEG velocity set and a velocity dispersion (index 
is); ii) the WIMP mass (labeled as mw) and 0; iii) the galactic halo model and all 
the other parameters involved in the calculation (index im)-, the likelihood ratio as a 
function of psgr can be defined: 



'{Ps 



"^aa;^(s,,so) [L'^;^'iPsgr = 0)] 



'''^(SI.SD) 



'(a 



sgr ) 



(10) 



where maxa.si ^o) l-^mC' (Psgr)] is the value of the likelihood function maximized with 
the respect to the particle cross sections. 

The functions Y.^^'{psgr) = —2ln{Xl^^') are asymptotically distributed as a 
chisquare with 1 degree of freedom. Some examples are given in Fig. 7 where three 
Y functions are plotted for some halo models and for some particle masses in the par- 
ticular case of a pure SI candidate and a SagDEG stream with velocity set Vsph and 




COl Q.C2 O.OS 0-04 0,05 0,06 0,07 0,08 0,09 

p„, (GeV/cm') 



Figure 7: Effect of the SagDEG contribution in the data fitting for 3 illustrative models. 
The example uses: SI candidate, a SagDEG stream with velocity set Vsph and velocity 
dispersion fo,sgr = 40 km/s; all the parameters for form factors and for the quenching 
factors are fixed at case A of ref. [3] (see also sec. 3.1 in the text). The considered 
halo models are: i) NEW halo (a = 1, /? = 3, 7 = 1, a = 20 kpc, A5 of [3]), vq = 220 
km/s, phaio = 0.74 GeV cm~^ and mw = 10 GeV (dotted line); ii) Evans' logarithmic 
halo [Re ^ kpc, q = l/\/2, CI of [3]), vq = 170 km/s, phaio = 0.56 GeV cm^^ and 
mw = 22 GeV (solid line); iii) Evans' logarithmic counter-rotating halo {Re — 5 kpc, 



q = I/V2, C2 of [3]), 
GeV (dashed line). 



Vq = 170 km/s, Phaio = 0.67 GeV cm , ij = 0.64 and mw — 20 
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velocity dispersion wo,sgr — 40 km/s. 

In particular, this figure shows three representative cases: i) a model where the 
SagDEG contribution worsens the data fit (dotted line); ii) a model where the SagDEG 
contribution improves the data fit providing a C.L. better than 3 a (solid line); iii) a 
model where the SagDEG contribution improves the data fit providing a C.L. lower 
than 3 a (dashed line). 

In order to investigate the presence of SagDEG, in all the considered halo models 
and adopted parameter uncertainties, for simplicity here we alternatively investigate 
only the purely SI and the purely SD cases, respectively. In the following - for each 
considered mw and is - the 90% C.L. allowed intervals on psgr are constructed requir- 
ing that [39]: 

>;^/= iPsgr) < rmn^^^^p^^^^ {Y^^-^^ (p,,,)) + 2.71. (11) 

In Fig. 8 the SagDEG density psgr allowed at 90% C.L. for pure SI coupling is shown 
as a function of mw- In Fig. 9 the same is shown for a SD coupled candidate in the 
particular cases of 6* = 0,7r/2,7r/4 and 2.435 (pure Zq coupling). 

From the figures 8 and 9 upper limits on the SagDEG density can be inferred. In 
particular, for some WIMP masses and for some halo models, these limits are compa- 
rable or improve the limit already given in ref. [20] (of the order of 0.07 GeV cm~'^) 
on the basis of considerations on AI/L. 

Moreover, figures 8 and 9 suggest that intervals not including psgr = at 90% C.L. 
exist for some values of mw- This points out a slightly preference for the presence 
of a SagDEG contribution in the data. However, considering the uncertainties on 
the SagDEG velocity and velocity dispersion (that is e.g. superimposing the allowed 
regions in the figures 8 and 9) in most of the considered scenarios the absence of 
SagDEG is still allowed at 90% C.L.. 
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Figure 8: SagDEG density psgr allowed at 90% C.L. (hatched area) for pure SI coupling 
as function of m^i/ values. Left panel: case of Vg* velocity set. Right panel: caseofV^sp/i 
(descending hatched) and Voti (ascending hatched) superimposed. The used stream 
velocity dispersion, wo,sgr values are: 20 km/s (dashed), 40 km/s (solid) and 60 km/s 
(dotted). 
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Figure 9: SagDEG density psgr allowed at 90% C.L. (hatched area) for pure SD 
coupling as function of mw values for some of the possible 9 values. Upper 4-box: 
case of Vs* velocity set. Lower 4-box: case of Vsph (descending hatched) and Vou 
(ascending hatched) superimposed. The used stream velocity dispersion, uo,sgr values 
are: 20 km/s (dashed), 40 km/s (solid) and 60 km/s (dotted). 
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It is worth to note that in many of the analyzed configurations the inclusion of 
the SagDEG contribution improves the data fit. For example, for SI candidate in the 
case of a stream with velocity set Vsph and velocity dispersion vo^sgr = 40 km/s, about 
67% of the configurations have an improvement of the data fit by the inclusion of the 
SagDEG; in particular, the improvement of about 18% of them is better than 2 a. 
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Figure 10: Example of the cumulative percentage distribution of psgr best-fit values 
providing a C.L. better than 2 a with the respect to the absence of SagDEG. A pure 
SI candidate and fixed SagDEG stream with velocity set Vsph and velocity dispersion 
VQ^sgr = 40 km/s have been considered. About 60% of these psgr best-fit values are 
below 0.1 GeV cm~'^. See text for implications. 



Another interesting information can be inferred by studying the pagr best-fit val- 
ues achieved for the various considered models. For this purpose, the cumulative 
percentage distribution of psgr best-fit values providing a C.L. better than 2 a with 
the respect to the absence of SagDEG is shown in Fig. 10. A pure SI candidate 



and fixed SagDEG stream with velocity set Vsph and velocity dispersion vq. 



sgr 



40 



km/s have been considered here as an example. About 60% of these models gives 
Psgr best-fit values below 0.1 GeV cm~^; in addition, the distribution peaks around 
psgr ~ 0.04 GeV cm~^. These latter values are intriguing considering the expectations 
on the stream density at Sun position - that is few % of the local dark halo - based 
on some theoretical studies about the disruption of the satellite galaxies falling in the 
Milky Way halo [40]. 

This preliminary analysis offers hints on the possibility to investigate halo fea- 
tures by annual modulation signature already at the level of sensitivity provided by 
DAMA/Nal. 
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5 Conclusion 

In this paper a preliminary study on the effect of the presence of Dark Matter particle 
streams in the galactic halo has been analysed on the basis of the annual modulation 
data collected by DAMA/Nal. 

In particular, the case of the Sagittarius Dwarf Elliptical Galaxy (which presently 
is the better known case) has been discussed here showing its effect on the allowed 
volumes/regions for some astrophysical, nuclear and particle physics scenarios related 
to the case of WIMP candidates. 

The potentiality of a similar approach to investigate the halo composition has 
also been pointed out as well as the possibility to derive experimental bounds on the 
possible contribution of the SagDEG to the local dark matter density. For some of the 
investigated WIMP masses, the order of magnitude of these bounds, obtained by local 
measurements ^, is in agreement with the existing bounds based on non-local M/L 
ratio observations. 

Other candidates and other non thermalizcd substructures will be addressed in near 
future studies; in particular, the availability of larger exposures by DAMA/LIBRA will 
offer the possibility of more efhcient discrimination capability among different possible 
scenarios. 
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